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Анализ свойств наноповерхностей 
по АСМ изображениям 


В данной статье представлена схема алгоритма анализа наноповерхности по изображениям атомно-силовой 
микроскопии. В основе анализа лежит коррекция неравномерного контраста и выделение объектов с округлой 
формой. Данный алгоритм позволил определить не только классические характеристики описания нано- 
поверхности, такие как шероховатость, но и дополнительные характеристики нанообъектов, позволяющие 
получить как их индивидуальное, так и групповое описание. 


Введение 


Атомно-силовая микроскопия позволяет получать достоверную информацию о то- 
пографии рельефа объекта [1]. При исследовании характеристик нанообъектов рас- 
сматриваются не только отдельные атомы или молекулы, но и их достаточно крупные 
макрообъединения. 

Для дальнейшей обработки таких изображений необходимо выделение объектов и их 
характеристик. Эта задача является довольно сложной, так как наноизображения имеют 
много особенностей. Ранее выполнялись попытки обработки подобных изображений [2-5], 
однако в данном случае каждый конкретный случай требует отдельного рассмотрения, так 
нами ранее обсуждалась обработка магнито-оптических изображений [6], [7]. В данной 
статье производится дальнейший анализ магнито-оптических изображений на примере 
изображений нанообъектов. 

В данной работе рассматривается использование универсальных программных па- 
кетов инженерного анализа — САЕ-систем при подготовке данных для проведения чи- 
сленных экспериментов по определению характеристик контактного взаимодействия 
нанообразцов и приводятся результаты вычислительных экспериментов по определению 
характеристик контактного взаимодействия нанообразцов. В работе приводятся результаты 
вычислительных экспериментов по определению характеристик их контактного взаимо- 
действия средствами пакета Г/5-РУМА. 


Особенности наноповерхности изображения 


В поставленной задаче определения шероховатости объекты представлены полу- 
сферами. На полутоновых наноскопических изображениях они представлены округлыми 
регионами с увеличивающейся яркостью к центру (рис. 1). 
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Рисунок 1 — Полутоновое изображение наноповерхности 
со сферическими объектами 
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Исходными данными являются полутоновые изображения, имеющие ряд дефектов. 
Горизонтальный профиль яркости исходного изображения представлен гладкой 
линией (рис. 2). 
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Рисунок 2 — Горизонтальный профиль яркости наноскопической поверхности 
по горизонтали 


По вертикали он имеет профиль с большим контрастом и небольшими выбро- 
сами (рис. 3). 
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Рисунок 3 — Вертикальный профиль яркости 


В результате на изображении присутствует неоднородный геометрический контраст, 
где составляющая изменений по вертикали значительно превышает составляющую по 
горизонтали. Этот эффект связан с условиями формирования изображения. Процесс по- 
строения изображения основан на сканировании поверхности кантилевером. В общем 
случае это позволяет получить трёхмерное топографическое изображение поверхности с 
высоким разрешением, которое представлено в виде полутонового изображения. Сканиро- 
вание выполняется по линиям. В результате, время, затраченное на получение соседних 
пикселей по вертикали, значительно превышает время получения пикселей по горизонтали. 
За это время на результат «успевает» повлиять целый ряд внешних помех: геометрии- 
ческий дрейф — поверхность смещается и сканирование выполняется не по соседней 
области, дрейф яркости — механический способ движения иглы чувствителен к любым 
колебаниям в окружающем пространстве. Кроме того, при горизонтальном сканировании 
для кантилевера свойственна некоторая инерционность. Это приводит к тому, что свойст- 
ва яркости каждой линии получают уникальные свойства. Такая неоднородность отра- 
жается на глобальной гистограмме яркости (рис. 4). Гистограмма наноповерхности, не- 
смотря на гауссоподобную форму, имеет много отдельных импульсов, которые резко 
отличаются от своих соседей. Все эти дефекты влияют на качество обработки изобра- 
жения и требуют их удаления или компенсации. 


Рисунок 4 — Гистограмма изображения наноповерхности 
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Коррекция локальной неоднородности контраста 


Геометрическая неоднородность контраста на изображениях, получаемых методами 
атомно-силовой микроскопии, является серьезной проблемой, решение которой в общем 
случае до сих пор не найдено. В данной работе мы предлагаем два подхода: первый за- 
ключается в генерации возможного фона для изображения, второй основан на вычислении 
характеристик локального контраста для вертикальных и горизонтальных линий. Оба под- 
хода основываются на предположении о симметричности объектов на наноповерхности. 

В первом способе генерируется два фоновых образа. Один создается много- 
кратным применением простой усредняющей свертки с большим растром (для 
изображения 256 х 256 - растр 14 х 14): 
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где /.,— значение пикселя на изображении с координатами х, у, п — полуразмер 
растра свертки. 
В результате такого преобразования получается изображение, соответствующее 


темному фону (рис. 5 6). ПЕНИЕ | 
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Рисунок 5 — Изображение наноповерхности (а) 
и сгенерированный для него темный фон (6) 


Второй фоновый образ соответствует максимальному уровню, т.е. «засветке 
изображения». 

Данное изображение генерируется несколькими итерациями полутоновой дилатации 
(для изображения 256 х 256 — 2 итерации округлым растром 5 Х 5). Дилатация множества А 
по множеству В (или относительно В) определяется как 

д®в= {| (В), ПА=8}. 

В основе этого соотношения лежит операция получения центрального отражения 
множества В относительно его начала координат (которое для краткости будем называть 
центром В) и затем сдвиг полученного множества в точку 2. При этом дилатация мно- 
жества А по В — это множество всех таких смещений, при которых множества В и А 
совпадают. Для полутонового пикселя происходит смещение границ, т.е. текущему 
пикселю присваивается значение его соседей. 

В результате получается изображение слившихся объектов (рис. 6). 


я 
Рисунок 6 — Изображение фона максимальной яркости для наноповерхности 


248 «Искусственный интеллект» 3’2011 


Анализ свойств наноповерхностей по АСМ изображениям 


Имея два предела для изображения по яркости, можно приблизительно определить 
локальный контраст каждой точки. Он приблизительно будет равен разнице между ними. 
На основе этих данных можно определить значения изменения контраста в каждой точке, 
но данная операция не является эффективной. Основываясь на наличии двух фоновых 
изображений, можно выполнить операцию коррекции фона, которая заодно компенси- 
рует недостатки пространственной неоднородности контраста: 


С _ Ч» В») — В,,) 
и В) 


где /.,-— значение яркости пикселя изображения с магнитооптической компонентой; 
В,,-— значение пикселя для изображения с постоянной компонентой; И’,, — значение 
пикселя для изображения при полном насыщении (засветке); И’. — максимальное зна- 
чение яркости изображения при насыщении яркости; С..,-— новое значение яркости пикселя 
в скорректированном изображении (рис. 7а). 
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Рисунок 7 — Изображение наноповерхности после коррекции фона (а), 
после дополнительной медианной фильтрации (б) 

и после локальной коррекции контраста (в) 

Однако на полученном изображении остались горизонтальные линии, и их контраст 
даже увеличился. Подобный недостаток компенсируется медианной фильтрацией, позво- 
ляющей сгладить поверхность изображения. 

Второй способ заключается в вычислении локального контраста. Его суть состоит в 
определении числового значения локального контраста для каждой линии изображения, 
нелинейном его усилении и восстановлении той же линии изображения с измененной 
яркостью, что обеспечивает в сравнении с исходным изображением усиление локального 
контраста по линиям. Структурно процедура усиления локального контраста состоит из 
трех основных этапов и используется для каждого элемента [(1/) с координатами (1) 
исходного изображения Г. 

На практике эффективность способа недостаточна для обработки изображений, ко- 
торые содержат мелкие детали — зазоры между объектами. Причина заключается в том, 
что локальный контраст определяется по формуле, где его значение пропорционально 
мере отличия центрального элемента изображения от окружающего фона по значению 
яркости. Составными элементами этой формулы являются непосредственные значения 
элементов или их усредненные значения, что приводит к неполному описанию текстуры 
локальной области. Наиболее полно такие характеристики текстуры как однородность, 
«шершавость» и зернистость описываются статистическими методами. 

Одним из наиболее простых методов описания текстуры является использование 
моментов гистограммы интенсивностей элементов изображения. Пусть [, — случайная 
величина, которая определяет дискретную интенсивность изображения, Н([(1))) — со- 
ответствующие значения гистограммы. Известно, что и-й момент /(1/) относительно 
среднего значения определяется формулой 


и,(р)= Хе.) Н,)), 


где Г -— среднее значение яркостей элементов локальной окрестности Й'. 


«Штучний 1нтелект» 32011 249 


Недзьведь А.М., Абламейко С.В., Недьзведь О.В. 


Из этого выражения следует, что до = 1, а и! = 0. Второй момент соответствует 
дисперсии о”([) и служит для описания текстуры. Он является также мерой контраста 
интенсивности и применяется для описания однородности поверхностей. В результате в 
качестве меры контраста текстуры предложено использовать выражение 


И = 7 , 


где о”([)— дисперсия в окрестности и х т, к = 0,8 — коэффициент нормирования. С(/) 
равно нулю для окрестностей с постоянной интенсивностью и единице — для больших 
значений о”(Г), что полностью отвечает требованиям определения локального контраста. 

Таким образом, в предложенном методе на его первом этапе для каждого элемента 
изображения необходимо вычислить локальный контраст как С(1.). Далее осуществляется 
нелинейное преобразование локального контраста С(1.]), и затем восстанавливается изобра- 
жение путем определения нового значения яркости пикселей линии [. (1,]), для чего исполь- 


зуется выражение: 
0,5 
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Описанная процедура повторяется для каждого элемента каждой горизонтальной 
линии изображения. 

В предложенном алгоритме используется статистическое определение локальных 
контрастов, благодаря чему учитываются как однородность текстуры, так и ее «Шер- 
шавость» и зернистость (рис. 7 в). 

Эксперименты показали, что результат, полученный вторым способом, является 
более качественным, в смысле своей четкости, однако первый способ работает в три раза 
быстрее. 


ЕЛ =ЫЬЛ+ 


Сегментация сферических объектов на наноповерхности 


Основой для определения характеристик поверхности (размера, формы, высоты и 
положения) являются образующие ее объекты. Для определения этих характеристик надо 
выделить соответствующие объекты на изображении. 

Выделение объектов может осуществляться разными способами. Нами были разра- 
ботаны два способа выделения объектов. Исходя из предположения о том, что объекты 
представлены сферами и пересекаются, их можно представить как топологическую по- 
верхность. В этом случае можно определить положение впадин между ними. 

Для определения ограничивающих линий между объектами существует алгоритм 
водораздела. Этот алгоритм является основным инструментом математической морфо- 
логии для сегментации изображений. Однако полученное изображение часто получается 
чрезмерно сегментированным на большое число маленьких областей (бассейнов), боль- 
шинство из которых не является значимыми при решении проблемы. 

В топографическом представлении изображения, численные значения (например, 
уровни серого) каждого пикселя выступают в качестве высоты этой точки. Преобразование 
водораздела вычисляет водосборные бассейны и линии хребтов. При этом в качестве 
водосборных бассейнов рассматриваются соответствующие области изображения, а 
линии хребтов являются границами этих областей. Основной проблемой данного алгорит- 
ма является чрезмерная сегментация, поскольку все границы и шумы отображаются в 
градиенте, что делает необходимым процесс удаления шума. В результате прямого при- 
менения этого алгоритма к изображению наноповерхности получаются дополнительные 
ложные объекты (рис. 8). 
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Рисунок 8 — Результат работы алгоритма водораздела: а) примененный 
непосредственно к изображению наноповерхности, 6) примененный к изображению 
наноповерхности после фильтрации шумов, в) примененный к изображению 
наноповерхности после коррекции локального контраста и фильтрации шумов 


Улучшение результатов работы алгоритма водораздела заключается в сглаживании 
поверхности изображения. Для этого обычно применяются различные фильтры: линей- 
ный усредняющий, нелинейный медианный и морфологический. Однако в случае нано- 
поверхности их применение перед преобразованием водораздела приводит к потере 
объектов на изображении, так как некоторые объекты сливаются вместе. 

Такой результат получается из-за геометрической неоднородности контраста, о ко- 
торой говорилось в предыдущей части. Применение коррекции контраста перед фильтра- 
цией, удаляющей шумы, значительно улучшает результат (рис. 8 в). 

Алгоритм позволяет получить границы объектов, но линии раздела получаются 
немного рубленые, так как преобразовании используются городские метрики вычисления 
расстояний. 

Как говорилось выше, объекты на изображении представляют собой полусферы, 
однако изображение границ преимущественно состоит из многоугольников. Для решения 
этой проблемы разработан специальный алгоритм коррекции формы на основе математи- 
ческой морфологии. 

Первым шагом алгоритма является инверсия границ (рис. 9 а). 


Рисунок 9 — Процесс бинарной обработки а) инвертированные границы объектов, 
6) результат удаления объектов прижатых к краю, в) удаление мелких объектов, 
г) коррекция формы объектов 


Затем удаляются объекты, прижатые к краю (рис. 9 6). 

На изображении присутствуют мелкие объекты, которые возникли на пустых участ- 
ках. Они удаляются на основе фильтрации по площади (рис. 9 в). 

Затем выполняется коррекция границ объектов с помощью операции морфологи- 
ческого отмыкания (рис. 9 г). 

Данный алгоритм имеет ряд преимуществ по сравнению с первоначальным спосо- 
бом получения радиусов, так как позволяет вычислить ряд геометрических и яркостных 
характеристик, которые могут использоваться для более полного описания. 
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Вычисление характеристик объектов 


Описание поверхности может быть расширено за счет использования обобщен- 
ных характеристик объектов. Используя такие характеристики можно определить как 
свойства индивидуальных объектов на поверхности материала, так и свойства всей по- 
верхности. 

Объекты определяются как области на изображении, соответствующие информа- 
тивным структурам. Для извлечения полной информации об объектах изображения и исполь- 
зования ее для постановки диагноза необходимо вычислить характеристики (параметры) 
объектов изображения. На их основе делаются выводы о типе объекта или о качестве со- 
держимого изображения. В зависимости от объектов определяют способ его описания и 
его характеристики. Основная классификация объектов определяется по их размерности. 
В области анализа изображений наиболее распространены три размерности, однако в 
частных случаях может использоваться и больше. В зависимости от размерности опре- 
деляются классы характеристик, которые описывают объект. 

Работа с трёхмерными объектами подразумевает три непременных этапа: построение 
изображения (визуализация), получение описательных характеристик, хранение инфор- 
мации об изображении. При этом первых два этапа в значительной степени взаимозави- 
симы. Процесс визуализации во многом связан с предварительным анализом основных 
характеристик объекта. Процесс визуализации и форматы сохранения полученных данных 
были подробно описаны на предыдущем этапе выполнения задания. Сейчас наша задача 
состоит в описании параметров, характеризующих объёмные особенности поверхности 
микроструктур, а также методов, использованных для получения этих характеристик. 


Выводы 


В данной статье представлена схема алгоритма анализа наноповерхности по изобра- 
жениям атомно-силовой микроскопии. В основе анализа лежит коррекция неравномер- 
ного контраста и выделение объектов с округлой формой. Данный алгоритм позволил 
определить не только классические характеристики описания наноповерхности, такие как 
шероховатость, он и дополнительные характеристики нанообъектов, позволяющие полу- 
чить как их индивидуальное, так и групповое описание. Данный алгоритм реализован в 
программном комплексе ТР (рис. 10). Результаты вычислительных экспериментов для 
определения характеристик контактного взаимодействия были проведены на основе 
анализа реальных изображений металлов и керамики. Полученные результаты покА- 
зывают, как на изображениях наноповерхности для объектов можно вычислить ряд 
характеристик и особенностей даже в таких сложных изображениях, как изображения 
наноструктур. Дополнительные характеристики позволили описать новые физические 
свойства наноповерхности без дополнительных дорогостоящих испытаний. 


018. 
3.4641, 2.73861, 1.26491, 1.33333, 


Рисунок 10 — Копия экрана программы ТПР для анализа наноповерхности 
Работа частично сопровождалась проектами МНТЦ В-1682 и В-1636. 
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А.М. Недзьведь, С.В. Абламейко, О.В.Недзьведь 

Анал!з властивостей наноповерхностей за АСМ зображенням 

У дашй статт! представлена схема алгоритму анал1зу наноповерхност! за зображеннями атомно-силово! 
микроскопй. В основ! аналзу лежить корекшя нервномрного контрасту 1 видлення об’екив з округлою 
формою. Даний алгоритм дозволив визначити не ‘ильки класичн! характеристики опису наноповерхност, 
так! як шорстксть, але 1 додатков1 характеристики нанооб’екив, що дозволяють отримати як 1х 1ндивдуаль- 
ний, так 1 груповий опис. 


А.М. Меасуеа, 5.7. АМатеуко, О.Т. №е4зуе4а 

Апа1у$15 оГ Мапо-Зиг асе Ееафиге$ Бу АЕМ Ппаде Ргосез$т? 

ш 15 рарег Фе зсВеште оЁ а!еогИ бт оЁ папо-зигРасе ипазе апа[уз$1$ оЁ аюпис Ктгсе плсгозсору \а5 
Чезсге4. ТБе апа[у$1$ 1$ Базе оп соггесНоп о{ ипеуеп сопгаз( ап обес з@есНоп \\/Ий а гоипае4 Варе. 
ТЬ$ а!оогбт аПо\уе чз ю деепише поё ошу е с1аз$1са| свагасетзИс$ оЁ папо-зитРасе зисН а$ 
тоиеВпезз, Бай и шсш4е4 а опа] сБагас{ет15Нс$ оЁ зеуега| папо-об]есв. 


Статья поступила в редакцию 22.06.2011. 


«Штучний 1нтелект» 32011 253 


